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Abstract 

Measurement-based quantum computing (MBQC) is a model of quantum computing 
that proceeds by sequential measurements of individual spins in an entangled resource state. 
However, it remains a challenge to produce efficiently such resource states. Would it be 
possible to generate these states by simply cooling a quantum many-body system to its 
ground state? Cluster states, the canonical resource states for MBQC, do not occur naturally 
as unique ground states of physical systems. This inherent hurdle has led to a significant 
effort to identify alternative resource states that appear as ground states in spin lattices. 
Recently, some interesting candidates have been identified with various valence-bond-solid 
(VBS) states. In this review, we provide a pedagogical introduction to recent progress 
regarding MBQC with VBS states as possible resource states. This study has led to an 
interesting interdisciplinary research area at the interface of quantum information science 
and condensed matter physics. 



1 Introduction 

Although Richard Feynman ^ mooted the idea of a quantum computer nearly thirty years ago 
in 1982, serious research work into quantum computing probably only appeared about fifteen 
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years ago when Peter Shor demonstrated that a quantum computer, if available, could perform 
prime number factorization exponentially faster than a conventional computer which we typi- 
cally call "classical computer". A quantum computer essentially exploits inherently quantum 
mechanical properties like the linear superposition principle and entanglement to improve com- 
putational speed and complexity. 

In Feynman's original idea, a quantum computer, if constructed, would be also an ideal 
machine for simulating naturally-occurring many -body quantum problems. After all, nature is 
largely quantum mechanical and the ideal platform for solving problems in nature should be 
a quantum mechanical device. From a theoretical computer science perspective, the rigorous 
foundation of a quantum computing device was laid, in 1985, by David Deutsch at Oxford who 
cleverly formulated a theoretical model called the quantum Turing machine [i2J, which has been 
a fundamental lego block for the study of a quantum computer. 

In 1989, Deutsch introduced another model of fundamental importance for quantum com- 
putation: the quantum circuit model [3]. Yao then showed in 1993 that the circuit model is 
actually equivalent to the quantum Turing machine in terms of their computational power 
As we have mentioned before, Shor published his famous algorithm for efficient prime numbers 
factorization on a quantum computer in 1994 (Si. Factoring is believed to be a difficult problem 
to solve on a classical computer and it lies at the heart of modern cryptographic technology: 
the Rivest-Shamir-Adleman (RSA) scheme. Shor's theoretical insight, if realized, thus poses 
an immediate threat to classical cryptographic security. Equipped with quantum computers, all 
security systems currently in use based on the computational difficulty of large prime number 
factorization would be rendered useless. 

The basic entity in classical computing is the bit. A bit is usually encoded in several ways, 
for instance, the voltage or current in a wire or the energy in a capacitor. In a quantum com- 
puter, the information can also be encoded in different ways: for instance, the polarization of a 
photon or the energy levels of a two-level atomic or molecular system. Unlike classical bit, the 
information for a two-level quantum system can be encoded in a superposition of the two levels 
of the system. Such encoding is called a quantum bit or a qubit. Thus, qubits form the basic 
building blocks for a quantum computer. 

Classical computation proceeds essentially through a series of wires and gates. The quantum 
circuit model in which algorithms are implemented through a series of quantum gates thus 
remains the most promising paradigm for the large scale realization of a quantum computer. 
For a quantum computer, the circuit model applies a sequence of unitary operations on an 
initial quantum state of a system with n qubits and evolves the quantum system according 
to the Schroedinger equation. Although each single gate usually acts non-trivially on only one 
or two qubits, the concatenation of these gates proves to be sufficient for producing any unitary 
operation U on the n-qubit initial state a property known as "universal" computation. The 
application of this sequence of gates transforms an initial state \'ipi) into a final state \ipf) = 
U\%jji). The final computational result is then extracted from appropriate measurements on the 
final state. 

The realization of the quantum circuit model in various physical systems has been exten- 
sively studied. These approaches include implementation through nuclear magnetic resonance 
H, the manipulations of atoms or ions in ion traps Q, the manipulation of neutral atom 
(HI, implementation with cavity QED f9], the optical platform with linear or nonlinear optical 
devices IfTOll . manipulation of electrons or atoms in the solid state [fTTTl . superconducting de- 
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vices [[T2II . and a "unique" qubit scheme [|T3ll . Regardless of the approach in the realization of 
large scale quantum computing, DiVincenzo has elegantly summarized five important criteria 
for the physical implementation of future quantum computers |[T4ll . and each one of these cri- 
teria has been carefully examined (71 [H |9l \Wi [HI [I2l [IS [TH. To date none of these systems 
is capable of realizing a large scale quantum computer in the foreseeable future without any 
glitches: each system presents uniquely its own set of challenges and problems. And indeed one 
important stumbling block is that many quantum gates entangling two qubits cannot be imple- 
mented with high fidelity in many of these systems ifTSl . To overcome some of these limitations, 
researchers have sought other paradigms of quantum computation. Although these alternative 
paradigms are equivalent to the quantum circuit model in terms of computational power, they 
are very different in terms of real physical realization. Among these other paradigms, some of 
more promising ones are measurement-based quantum computing [|T6ll . topological quantum 
computing ifTTl and adiabatic quantum computing lITSl . 

In this review, we focus on the measurement-based quantum computing (MBQC) model, 
first introduced by Raussendorf and Briegel in 2001 lfT6]| . MBQC is a model of quantum com- 
puting that is equivalent to the circuit model in terms of computational power, given that it can 
simulate every single-qubit gate and two-qubit gate in the circuit model in an efficient way, and 
it is thus universal for quantum computing. However, the physical realization of MBQC is very 
different from that of the circuit model. It starts from a highly entangled multi-qubit state \^c), 
called a cluster state. Sequential single-qubit measurements are then performed on the cluster 
state to realize quantum computing. The advantage of MBQC is that once I'ifc) is prepared, 
no entangling gates are needed in process of the computation. Or one may imagine that all the 
entanglement needed for a quantum computation is already embodied in |^c)- MBQC thus 
seems very appealing for the implementation of quantum computing provided that can be 
efficiently prepared initially and maintained during the process of measurements. 

Perhaps the most natural and simplest means of preparing an entangled quantum state is 
to search for a many-body system in which the unique ground state is the required entangled 
resource. With correct cooling techniques, one should then be able to acquire the ground state 
and hence the entangled resource state needed for MBQC. Unfortunately, it turns out that |\l/c) 
does not naturally occur as unique ground states of any naturally-occurring physical systems. 
More precisely, it was shown by Nielsen that \^c) cannot even be a unique ground state of 
any Hamiltonian of qubits involving at most two-particle interactions |fT9ll . It is then highly 
desirable seeking for other many-body quantum states, with which MBQC can be performed 
in a similar manner as is performed on I'^c), where only single-particle measurements are 
needed to implement universal quantum computing. We call these alternative states for MBQC 
"resource states". Moreover, we want these resource states to naturally occur as unique ground 
states of some lattice spin systems, that is, unique ground states of Hamiltonians involving only 
two-particle nearest-neighbor interactions on a lattice. In addition, one also wants the resource 
states to be stable during the process of measurements, that is, to remain as a ground state of 
the Hamiltonian after the particles are measured and discarded. To satisfy this condition, it is 
sufficient that the Hamiltonian is "frustration-free", that is, the ground state of the Hamiltonian 
is also a ground state of each interaction term of the Hamiltonian. Another desired property 
is that the Hamiltonian is gapped, meaning that it has a constant energy gap regardless of the 
system size (i.e., the number of particles), so the ground state can be stable against thermal 
fluctuation. 
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In 2007, Gross and Eisert took the first step towards the construction of such realistic re- 
source states for MBQC EOl . They introduced an interesting framework for producing resource 
states which is closely related to the theoretical description of valence-bond-solids (VBS). The 
gist of their proposal is that even if the quantum computation is actually performed on a system 
of qubits, the actual physical system is not just a system of particles with information encoded 
in two different levels. These particles could ideally possess higher spins - a feature which 
greatly facilitates the search for realistic lattice Hamiltonians. Their framework allows one to 
construct various resource states with different spin values and on different types of lattices. 
For resource states associated with realistic Hamiltonians, they found Hamiltonians on a spin-1 
chain, involving only nearest-neighbor interactions, whose unique ground states are resource 
states for MBQC. Moreover, these systems are also frustration-free and gapped. Unfortunately, 
their scheme applies specifically to one-dimensional(lD) spin chains. Since a ID spin chain can 
only process quantum computation of a single qubit, their framework are not universal. It was 
also not obvious how one could extend the scheme to find any realistic two-dimensional(2D) re- 
source states for MBQC. In 2008, Brennen and Miyake showed that a famous one-dimensional 
VBS state, the Affleck- Kennedy-Lieb-Tasaki (AKLT) state on a spin-1 chain IIITI with certain 
boundary condition, can act as a resource state for MBQC Il22ll . Although this is still a result 
for ID chain, it sharpens the connection between the studies of MBQC and VBS states [|23l . 

In 2009, Chen et al. took an important step towards the construction of 2D realistic resource 
states for MBQC [j24|. They showed that a state \ ^tric) of ^ spin-5/2 system on the honeycomb 
lattice, called the tri-Cluster state, is a unique ground state of a frustration-free, gapped, two- 
body Hamiltonian with only nearest-neighbor interactions. Although the form of the actual spin 
Hamiltonian appears complicated, their method does open up new possibilities for constructing 
other 2D realistic resource states for MBQC. In early 2010, Cai et al. then constructed a realistic 
resource state of a spin-3/2 system on the 2D octagonal lattice, called a quantum magnet, which 
is in turn based on proper coupling of spin-3/2 quasi-AKLT chains. In the same year, Chen et 
al. developed an interesting alternative viewpoint on the universality of VBS states for MBQC, 
by showing that these known resource states can be reduced to the cluster states via adaptive 
local measurements at a constant cost ll25l . This viewpoint turns out to be particularly useful for 
studying the universality of VBS states for MBQC, as cluster states themselves can be viewed 
as VBS states that was first realized by Verstraete and Cirac in 2004 [i26i1 . Also based on this 
viewpoint, Wei et al. reinterpreted the universality of the quantum magnet by reducing this 
resource states to a 2D cluster state ||27]| . Interestingly, Wei et al. and Miyake independently 
discovered that the 2D AKLT state on the honeycomb lattice can also be a resource state for 
MBQC |l28l|29l[30l|. Although it is still unclear if the 2D AKLT Hamiltonian on the honeycomb 
lattice is gapped OH . this result nicely links the study of MBQC and VBS. 

The "quick reduction" from spin-5/2 to spin-3/2 systems for realistic 2D resource states for 
MBQC naturally gave further impetus for the search of realistic resource states in a spin-1 /2 
system. Unfortunately, it was recently shown by Chen et al. that this may not be possible. 
Chen et al. provided a no-go theorem showing that there does not exist a unique ground state 
for a two-body frustration-free Hamiltonian that can simultaneously be a resource state for 
MBQC [[32ll for a qubit system. It was further shown by Ji et al. that indeed for two-level sys- 
tems, the structure of the ground-state space for any two-body frustration-free Hamiltonian can 
be fully characterized [33] and none of these states corresponds to a resource state for MBQC. 
These results herald bad news for the practical realization of MBQC based on VBS resource 
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states since spin- 1/2 systems appear to be the most prevalent systems in nature. One could try 
to overcome this situation by relaxing the requirement for "frustration-free" or "uniqueness" 
and resort to some other physical mechanisms, such as topological protection ll34l [35l or per- 
turbation Il36l l37l [381 . Finally it is worth noting that it is still an open question whether one 
can find a realistic spin-1 VBS resource state for MBQC with a two-body Hamiltonian which is 
both frustration-free and gapped. 

This review provides a pedagogical introduction to recent progress regarding MBQC with 
VBS states, which, we believe, is an interesting interdisciplinary research area at the interface of 
quantum information science and condensed matter physics. We assume no prior background 
on quantum computing of the readers, but we do assume basic knowledge of quantum physics, 
condensed matter physics and VBS states. For readers who are not familiar with VBS states, 
we would like to refer them to a recent review regarding entanglement in VBS states, which is 
published in the 2010 volume of this journal f23]. We will cover the following topics in this 
review. 

• A brief introduction to the circuit model of quantum computing. 

• An introduction to MBQC based on cluster states: what a cluster state is, how MBQC 
with cluster states simulates the quantum circuit model. 

• A viewpoint on the universality of the resource states for MBQC introduced in [|25ll . 

• MBQC in ID valence-bond chains, in particular, the universality of the ID AKLT chain 
for processing a single-qubit information. 

• The spin-5 / 2 tri-cluster state | "^tric) on the honeycomb lattice, its universality for MBQC 
and the properties of the corresponding tri-cluster Hamiltonian. 

• Resource states for MBQC in spin-3/2 systems, including the quantum magnets on the 
2D octagonal lattice and the 2D AKLT state on the honeycomb lattice. 

• The no-go theorem of resource states for MBQC in spin- 1/2 systems. 



2 The circuit model of quantum computing 

The basis for a two-level quantum system (qubit) is typically denoted by |0) and 1 1). A quantum 
operation on a qubit is a 2 x 2 unitary matrix, called a "single-qubit quantum gate". The Pauli 
matrices 

X=(° iV and Y=(° and Z^(i °V (1) 



together with the identity opertator 



.0 1 

form a basis for 2x2 matrices. The Hadamard gate H is given by 



/=(i ?) (2) 



V2 V 



H = ^{] \]. (3) 
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And other important single-qubit gates are the X, Y, Z rotations given by 



Xe = exp{-i9X/2), and Yg = exp{-ieY/2), and Zq = exp{-i9Z/2). (4) 

A basis for an n-qubit system is chosen as the tensor products of |0)s and | l)s. For instance, 
for n = 2, the four basis states are 

|0)® |0), |0)®|1), |0), |1), (5) 

which are in short written as 

|00), |01), |10), 111). (6) 

Unitary operations acting on two qubits are called "two-qubit quantum gates". The most- 
commonly used two-qubit quantum gate is the controUed-NOT gate, which takes \x) ® \y) to 
\x) \y (B x), where x,y E {0, 1} and © is the addition mod 2. Here the first qubit is called 
the control qubit, which remains unchanged, and the second qubit is called the target qubit, 
which is flipped if the control qubit is 1. In the basis of Eq.© the matrix of a controUed-NOT 
gate is then given by 



/I 0\ 

10 

1 

\0 1 0/ 



(7) 



Similarly, a controUed-NOT gate with the second qubit the control qubit takes |x) ® \y) to 

|x © (g) \y). 

Another important two-qubit gate is the controUed-Z gate, which transforms the basis of 
Eq.® in the following way: 

|00) ^ |00), |01) ^ |01), |10) ^ |10), 111) ^ -|11). (8) 

Given that the controUed-Z gate is symmetric between the two qubits, it is not necessary to 
specify which one is the control qubit and which one is the target qubit. 

In the circuit model of quantum computing, the initial n-qubit state \ipi) is usually chosen as 
the all |0) state |0)®|0) ■ ■ -©lO), which is in short written as |00 ■ ■ - 0) or |0)®". Then a sequence 
of single- and two-qubit quantum gates are applied on \ipi) to result in a final state [(pf). And 
finally single-qubit measurements are performed on each qubit, usually in the {|0), |1)} basis, 
to obtain the result of the computation. 

Indeed, if one is able to implement any two-qubit quantum gate on any two qubits, then 
any n-qubit unitary operation can be implemented by applying a sequence of two-qubit unitary 
gates |[39ll . In short, all two-qubit unitary gates together are universal for quantum comput- 
ing, meaning that they can be used to implement any quantum computation. Moreover, any 
two-qubit unitary gate can be implemented through a sequence of single-qubit unitary gates to- 
gether with the controUed-NOT (or controlled- Z) gate ||39]| . Therefore, single-qubit gates plus 
the controUed-NOT gate (single-qubit gates plus the controUed-Z gate) are also universal for 
quantum computing. 

To implement a particular algorithm, a circuit diagram is typically used to describe a quan- 
tum circuit. An example of a circuit diagram is given in Fig.[T] Each wire represents a qubit, and 
the time goes from left to right. In this example, there are total three qubits involved. The very 
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left side gives the initial state of the qubits, which is \ip) ® |0) ® |0), where \ip) is a single-qubit 
state. The very right side gives the final state of the qubits, which is Z\'mi) (g) H\m2) ® 
where mi,m2 E {0, 1} and \^') is some single-qubit state whose actual value depends on \ip) 
and 9. Each quantum gate in Fig.[T]is explained in Fig. [21 



1^) 
|0) 
|0) 



H 




Z\m,) 



(Hpi H\m,} 



Z, 



Figure 1 : A quantum circuit 



— U — 



(a) 



(&) 




Figure 2: (a) A single-qubit gate f/; (b) A controUed-NOT gate with the top qubit as the control 
qubit and the bottom qubit as the target qubit; (c) A controUed-Z gate; (d) A single-qubit mea- 
surement in the basis {f/|0), f/| 1)}, and the measurement result mi is obtained which will give 
an output state t/|mi). 



3 Measurement-based quantum computing 

This section provides an introduction to MBQC based on cluster states, first introduced by 
Raussendorf and Briegel in 2001 lfT6l . First let us briefly describe the way MBQC works. 
Fig. [3] provides a schematic illustration of MBQC. In this figure, each circle represents a qubit, 
which sits on a two-dimensional square lattice. Each pair of qubits that are linked by a solid 
line are called neighbors. The system of qubits is prepared to an initial cluster state |^c)- 
Then each qubit will be measured sequentially, from the left columns to the right columns, so 
the information flows from left to right. The arrow in each circle illustrates the direction of 
the spin that is measured, and the choice of directions are dependent on the results of earlier 
measurements. At the end of the procedure, all the spins are measured, and the measurement 
results together give the result of the computation. There is another name used for MBQC 
in literature, namely "one-way quantum computing", as the initial cluster state is completely 
destroyed after the computation. 

In Sec. 13. 1[ we will take a closer look at the cluster state. Then in Sec. 13.21 we show 
how MBQC can simulate any quantum circuit in an efficient way, hence universal quantum 
computing can be implemented by MBQC. 

3.1 The cluster state 

The cluster state was first introduced by Briegel and Raussendorf shortly before they introduced 
the model of MBQC HOll . The term "cluster state" actually refers to a family of quantum states. 
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information flow 



Figure 3: Measurement-Based Quantum Computing 

which are quantum states of qubits associated with certain lattices (or just simply arbitrary 
graphs). For any graph of n vertices, one can then define a cluster state. For examples, the 
graph can be the 2D square lattice discussed in Fig.[3l the 2D honeycomb lattice given in Fig.lH 
or a simple graph of three vertices given in Fig. [2 a). 

(><X><X><) 

Figure 4: The honeycomb lattice 

The creation of the cluster state associated with a given graph can be described, in terms of 
quantum circuit, in three steps. 

• Initialize every qubit in the state |0). 

• Apply a Hadamard gate on each of the qubits. 

• Apply a controlled- Z gate on each pair of qubits who are neighbors on the graph (i.e., 
whose corresponding graph vertices are connected by a solid line). 

As an example, a quantum circuit creates the cluster state associated with the graph of three 
vertices given in Fig. 13a) is given in Fig. Ob). 

For each vertex j in a given graph, denote the neighboring qubits of j by nb(j). Now 
we show that the cluster state \^c) associated with the graph is an eigenstate of the operator 
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(a) 



|0) 
|0) 
|0) 



H 

IT 
IT 

(b) 



Figure 5: (a) A three-qubit graph; (b) A circuit creates the cluster state associated with the 
three-qubit graph 



Xj <S)kenb(j) with eigenvalue 1, where Xj (Zk) is the Pauli X (Z) operator acting on the jth 
(kth) qubit. That is, 

Xj (g) Zfel^c) = l^c) (9) 

fcenb(j) 

We again take the simple example given in Fig. \5\ We denote the qubits from left to right 
as 1, 2, 3, and the controUed-Z gate acting on the jth and kth qubits by Sjk- Then the cluster 
state associated with the graph is given by (according to the circuit given by Ob)) 

\^c) = S2^S^2H^H,H,\000), (10) 

where Hj is the Hadamard operator acting on the j-th qubit. Now we show that Eq. ^ holds. 
That is, 



X,Z2\^c) 

To show this, observe the following identities: 



^c). (11) 



ZjSjk — SjkZj and XjSj^ — SjkXjZ^, 

ZjHj = HjXj and XjHj = HjZj. (12) 

Substituting Eq. (flOl) into Eq. (fT2l) will immediate result in Eq. (fTTI) . And it is straightforward to 
show that Eq. Q holds for the cluster state associated with any graph using a similar argument. 

One can see that the operator Xj (S)fcenb(j) has eigenvalues ±1, and the operators corre- 
sponding to different vertices commute. For instance, X1Z2 and Z1X2Z3 and commute 
with each other. Therefore, if we choose a Hamiltonian He as 

Hc = -Y,X, (g) Zfc, (13) 

3 fcenb(j) 

where the summation is running over all vertices of the graph, then | ^E'c) is obviously the unique 
ground state of the Hamiltonian He- Moreover, this Hamiltonian is obviously gapped (the 
entirely spectrum can be easily obtained) and frustration-free (as \'^c) is the ground state of 
each term in the summation). Note that such a Hamiltonian is in general not a two-body nearest- 
neighbor Hamiltonian. 
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3.2 Simulation of basic quantum gates 



In this subsection, we show that MBQC based on certain cluster states can simulate the quan- 
tum circuit model in an efficient way. As discussed before, to implement universal quantum 
computing, one only needs to realize single-qubit and controUed-NOT (or controUed-Z) gates. 
Apart from the original discussion on how MBQC simulates these gates lfT6l . there are other 
alternative discussions JHl |42l |43l IM IISI SB Here we would like to take the one pro- 
vided by Nielsen in |fT9l , which is in turn based on the circuit given in Fig.[6ta) that is proposed 
in 11471. 



1^) 
1+) 



m. 



(a) 



-X'^Hlij) 



1+) 



—X^HZe\^) 



Figure 6: (This figure is redrawn from Eqs.(ll) and (12) in |fT9l ). (a) A circuit for one-bit 
teleportation; (b) Generalized one-bit teleportation. 



The circuit in Fig.[6ta) is known as one-bit teleportation, where |+) = -^(|0) + To 
show how it works, \ci\ip) = a\Q) + . Then note 

i^i^i2(«|0) + /3|l))®-^(|0) + |l)) 

= -^(|0) ® H\,p) + XHm. (14) 

Measuring the first qubit in the {0, 1} basis gives the desired result. This circuit is then directly 
generalized to the one given in Fig. [6tb), as Zq commutes with Si2- 

Now we show that MBQC associated with the simple graph of three vertices given in 
Fig- Eta) simulates the circuit given in Fig. Hb). Here 1, 2 in each circle of the graph label 
the qubits 1 and 2. Then MBQC proceeds as follows: first measure qubit 1 in the basis of 
HZai\m), where m = {0, 1}. If m = is obtained, then measure qubit 2 in the basis of 
HZa^lm'), where m' = {0, 1}; if m = 1 is obtained, then measure qubit 2 in the basis of 
ifZ_a2|m'), where m' = {0, 1}. 




+) — HZa^ - HZa^— 



ib) 

Figure 7: (This figure is redrawn from Eqs.(13) and (14) in [fT9ll ). (a) A graph of three vertices; 
(b) The single-qubit quantum circuit that can be simulated by MBQC associated with the graph 
in (a). 



This procedure of MBQC can then be described by the quantum circuit given in Fig. Oa), 
which is equivalent to the one in Fig. (Hb). Now note here the circuit in each dashed box of 
Fig. (Mb) is nothing but the circuit of generalized one-bit teleportation given in Fig.[6tb). It then 
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1+)^ 

l+>- 
1+)- 




rn 



+ 



m 



Figure 8: (This figure is redrawn from Eqs.(13) and (14) in ifTgl ). (a) The circuit corresponds 
to the measurement procedure in Fig.[6ta); (b) The circuit equivalent to (a). 



follows that MBQC associated with the simple graph given in Fig. ITfa) simulates the circuit 
given in Fig.Hb). 

Note the ability to simulate the single-qubit gate of the form HZg for any 9 suffices to 
perform any single-qubit gate. This is because that HHZg = Zq and HZqH = Xq, so both 
X and Z rotations of any angle can be performed. And the idea for simulating the single-qubit 
gates by MBQC also generalizes to two-qubit gates. For example, MBQC associated with the 
graph given in Fig.^a) simulates the circuit of Fig.|9tb). The proof proceeds exactly along the 
same line as the description for MBQC associated with Fig.|3a), in which one can describe the 
procedure by a similar quantum circuits given in Fig. [8]^ a). 




HZa^ 






HZa^ 










HZ0, 






HZ0, 







ib) 

Figure 9: (This figure is redrawn from Eqs.(lO) and (17) in [[T9ll ). (a) A cluster state associated 
with a graph of six vertices; (b) The circuit corresponds to the measurement procedure in (a). 



The ability to simulate the single-qubit gate of the form HZq and the two-qubit gate given in 
Fig-Sb) thus can implement the controUed-Z gate. Together with arbitrary single-qubit gates, 
universal quantum computing can then be simulated by MBQC. 



4 Resource states for MBQC 

It has been shown in Sec. 13.21 that single-qubit gates and controUed-Z gates can be simulated 
by MBQC based on cluster states associated with certain kind of graphs. Universal quantum 
computing can then be implemented by putting these graphs together, which requires the graph 
to be of certain shape. Therefore not all cluster states can be used as resource states for universal 
MBQC, for instance cluster states associated with tree graphs (i.e., graph without cycles) cannot 
be resource states for MBQC [|48]| . However, many of them are possible, for instance, the cluster 
states associated with a square lattice or the honeycomb lattice [|49l . 
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The advantage of MBQC for the implementation of quantum computing in the real world 
is that once the resource state is prepared, then only single-qubit measurements are needed 
during the procedure of computing. It then remains a challenge how to efficiently produce such 
a resource state. To produce the cluster state, the direct way is to use the three-step quantum 
circuit discussed in Sec. 13. 1[ However, the third step requires applying controUed-Z gates 
on pairs of qubits. This then does not avoid the difficulty of implementing entangling gates, 
meaning one does not gain any advantage over the quantum circuit model. 

An appealing idea to overcome entangling operations is to seek physical systems in which 
one could obtain entangled resource states as unique ground states and the generation of these 
ground states is done through the cooling of quantum many-body systems. We already know 
that the cluster state is a unique ground state of the Hamiltonian He given in Eq. (|3]), and that 
He is gapped and frustration-free. However, He involve many-body interactions, not just two- 
body interactions and such Hamiltonians are generally not easy to find in practice. Also, the 
He associated with a 2D square lattice involves five-body interactions and the He associated 
with the 2D honeycomb lattice involves four-body interactions. Even the He associated with a 
ID chain involves at least three-body interactions. 

One therefore hopes that there exists a Hamiltonian with only two-body interactions that 
gives a particular cluster state as the unique ground state. However, this turns out to be impossi- 
ble as observed by Nielsen lfT9]| . The idea behind the proof is the following: if the Hamiltonian 
involves at most two-body interactions, then its ground-state energy is totally determined by the 
two-particle reduced density matrices of the ground states. Indeed, for any cluster state \^e) 
that is a resource state for MBQC, some of the eigenstates of He have the same two-particle 
reduced density matrices as l^c)- Therefore, if the cluster state |\l/c) is the ground state, then 
some other eigenstates of He are also ground states, meaning the cluster state cannot be the 
unique ground state of any two-body Hamiltonian. 

As a result, it is then highly desirable to seek resource states beyond the cluster states that 
are "naturally-occurring". We consider MBQC with a general spin-s system where s might be 
larger than 1/2, meaning we do not restrict ourselves to qubit or two-level systems. MBQC 
proceeds in a similar manner as in the qubit system, but now sequential single-particle mea- 
surements are performed on each spin-s particle, each then with 2s + 1 possible measurement 
outcomes. Of course, a spin- 1/2 system is highly desirable in practice, but simply by relaxing 
the restriction to spin-1/2 systems provides us some new ideas in the search for these highly 
entangled resource states. Accordingly, an "ideal state" l^j^) associated with a Hamiltonian Hid 
would satisfy the following conditions: 

1. lijjid) is universal for MBQC, i.e., \ipid) is a resource state for MBQC. 

2. lipid) is the unique ground state of Hid. 

3. Hid is a lattice spin Hamiltonian with only nearest-neighbor two-body interactions. 

4. Hid has a constant energy gap. 

5. Hid is frustration-free, meaning that \ipid) is the ground state of each interaction term of 

Condition 1—4 can be understood easily. Condition 5 ensures the stability of [tpid) throughout 
the process of MBQC, that is, after some particles are measured and discarded, the state of 



12 



the unmeasured particles continues to remain in the ground-state space of Hi^. Even if all five 
conditions are met, we would still like the value of spin s to be as small as possible. At a first 
glance, these demands appear daunting. Each condition alone is difficult to fulfill in general. 
However, we now argue that there are indeed simple examples which can reasonably satisfy all 
the conditions given. 

We begin with Condition 1. To discuss what kind of spin states are resource states for 
MBQC, one obvious way is to show that it can simulate the quantum circuit model in an efficient 
way. In practice this is not so easy to check. An important step was performed by Gross and 
Eisert in 2007 [20]. They introduced a framework for producing resource states, that is closely 
related to the theoretical description of valence-bond-solids (VBS) states and that allows one 
to construct various resource states with different spin values and on different type of lattices. 
Here in this review we take another viewpoint, which is also practical and readily applicable to 
VBS states. This main idea introduced by Chen et al. [(25l is the following observation: 

If a spin state \ip) can be reduced to a resource cluster state |\E'c) via adaptive local measure- 
ments at a constant cost, then is a resource state for MBQC. 

Moreover, this obviously sufficient condition for resource states might also be necessary. As 
shown by Chen et ah, all the known resource states do satisfy this condition. It turns out that 
this sufficient condition is convenient to apply for checking the universality for VBS states. To 
understand how this works, we would first mention another interesting observation: the cluster 
states can also be viewed as VBS states, which was first discussed by Verstraete and Cirac in 
2004 [|26ll. 

We start from a simple example of the cluster state associated with a ID chain graph, as 
shown in Fig. [lOl Here each circle represents a spin, called "physical spin." And each black 
dot also represents a spin, which we always choose as spin- 1/2 throughout this review, hence is 
called "virtual qubit." Each line connecting two virtual qubits is called bond. The state of the 
two virtual qubits connected by a bond is usually chosen as the singlet state which is 




Figure 10: A VBS state on a chain. 



In the usual description of VBS states, there is a projection Pj associated with the j-th 
physical spin, that projects the state of two virtual qubits inside the circle onto a d-dimensional 
subspace, hence the physical spin will be d-dimensional. For instance, for the famous Affleck- 
Kennedy-Lieb-Tasaki (AKLT) state ||2T]| . each Pj is chosen as the projection onto the symmetric 
subspace of the two virtual qubits, resulting physical a spin of dimension 3, i.e., spin-1. And 
the quantum state of the physical spins is then given by 

singlet )a, (16) 

j a 

where the subscript a labels the singlets. 
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We now show that the cluster state associated with a chain graph also possesses a VBS 
description, i.e., there exist projections P/s that project the state of two virtual qubits inside 
each circle onto a 2-dimensional subspace, such that the state given by Eq.( [T6l) is the ID 
cluster state. . For convenience, we choose the bond state as 

li^cbond) = ^(|00) + |01) + |10) - 111)) = H,Z,®X2\iPs^nglet), (17) 

which is equivalent to the singlet state by changing local basis of the virtual qubits, as shown 
by the second equality. 

Consider the projection is then chosen to be 

P=|0)(00| + |i)(ll|, (18) 

which is the same for each physical qubit. Here |0), |1) are the basis of the physical qubit. 
Recall that a cluster state associated with a line of n vertices is given by 

Note that 

/I \ 0" 1 

i+r^" =( -7^(10) + =^ E 1^1^- --n), (20) 

ii,i2,---,in 

where ii = {0, 1}. Therefore, for the term of Y[j 5'-,j+i|2ii2 • • • in), each string of the type 
I . . . 11 . . .) contributes a —1 phase factor. 

Now the VBS state of n physical qubits with the projection P given in Eq.ffTSl) as demon- 
strated in Fig.[TOl starting from n — 1 bonds given by Eq.([l7]), is given by 

\i^v) = l[Pj{\i^cl>ond)r''-'. (21) 

j 

To understand the resulted state, one can expand {\^cbond))'^"'~^ into a summation of terms of 
the form \i1i2 . . . «2n-2), with a phase factor either —1 or +1. After the projections Ylj Pj, only 
terms with the strings | . . . 00 . . .) and | . . . 11 . . .) will remain. Here the two qubits of the state 
00 or 11 correspond to any two virtual qubits in a same circle, which can then be replaced by 
and i after the projections as the states of the physical qubits. Furthermore, given the form of 
the bond state \iJcbond) as in Eq.(fT71). each string of the type | . . . 11 . . .) contributes a —1 phase 
factor. Therefore, the resulted state \'ipv) is nothing but the one-dimensional cluster state \^c)- 
Following a same argument, we choose the projections, instead, to be 

P= |0)(01| + |i)(10|, (22) 

which is the same for each physical qubit, then the resulted state \ipv) can also be the one- 
dimensional cluster state \^c) after performing some local transformation which maps |0) — )■ 
|i) and |i) |0). 

Let us now, consider, instead, the projection 

P' = |0)(00| + |i)(ll| + |2)(01| + |3)(10|, (23) 
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which is the same for each physical qubit. Note this projection actually projects each of the 
two virtual qubits in a same circle onto the full 4-dimensional Hilbert space. In other words, 
this projection does nothing but a relabeling of the four basis. Denote the resulted VBS state by 
iV'y), which is actually a state of a spin-3/2 system. 

Now consider the measurement on each of the spins of l^/'y), which is the projection onto 
either the {|0), or the {|2), |3)} subspace. Then according to the discussions above, the 
resulted state is a state of n qubits, which can be transformed to the ID cluster state |\l/c) after 
performing some local basis transformation dependent on each of the measurement results. In 
this way, Itp'y) reduces to a ID cluster state |\E'c) via adaptive local measurements at a constant 
cost, where one copy of [tp'y) will give one copy of \'^c)- Therefore, we know that \tp'y) is a 
resource state for MBQC of a single qubit. 

\ipy) thus provides an example of the sort of states described in [|25l regarding the usefulness 
of a quantum state for MBQC. The advantage of using \^y) instead of |\l/c) for MBQC of a 
single qubit is that \^y) can be associated with a Hamiltonian that satisfies all the Conditions 
2 — 5. This is obvious, as \4'y) is nothing but bunch of singlets by viewing each pair of the virtual 
qubits in a same circle as a spin-3/2 particle, whose corresponding Hamiltonian involves only 
two-body nearest-neighbor interactions, with \^y) the unique ground state, and is gapped and 
frustration-free. 

To generalize the above argument to 2D spin systems is just straightforward. For the cluster 
state associated with a 2D square lattice, as shown in Fig.dU one will just choose the projection 
of each four virtual qubits inside a same circle as 

P = |0)(0000| + (24) 

Similarly, the VBS state by simply viewing each four of the virtual qubits in a same circle 
as a 16-dimensional particle is a resource state for MBQC, which can be associated with a 
Hamiltonian that satisfies the Conditions 2 — 5. 




Figure 1 1 : A VBS state on a square lattice 



One can also consider the 2D honeycomb lattice. For the cluster state associated with the 
2D honeycomb lattice, as shown in Fig. [121 one will just choose the projection of each three 
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virtual qubits inside a same circle as 

P = |0)(000| + (25) 

Similarly, the VBS state by simply viewing each three of the virtual qubits in a same circle as 
a 8-dimensional particle (i.e., spin-7/2) is a resource state for MBQC, which can be associated 
with a Hamiltonian that satisfies Conditions 2 — 5. 




Figure 12: A VBS state on the honeycomb lattice 



We now have successfully constructed simple examples of "ideal states" which satisfy Con- 
ditions 1 — 5 as desired. The drawback of these constructions is obvious: each spin is of 
dimension 8 to perform universal MBQC. Even in the case for MBQC of a single qubit, one 
will need a spin-3/2 system. This is not practical for the realization of MBQC. So the most 
desirable thing is to reduce the spin value but still keep Conditions 1 — 5 satisfied. Compared to 
the simple construction of "ideal states" provided in this section, states with smaller spins turn 
out to be more difficult to construct, which requires many more other techniques apart from the 
sufficient condition for universality discussed above. This will be the topic for the rest of this 
review. We will start from the ID situation in the next section. 



5 MBQC in ID valence-bond chains 

In this section we consider ID VBS states as described in Fig.[TOl For convenience we would 
like to choose the bond state to be 

= ^(|00) + 111)) = ^1 ® X2\iJsinglet), (26) 

which is equivalent to Ifpsingiet) and \ipcbond) by changing local basis. The advantage of writing 
the bond state this way is the following. For each spin of the system corresponding to two 
virtual qubits inside the circle, let us write the projection as 

P = J2\'^)J2Mm]{kll (27) 

m k,l 
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where k,l = {0, 1}. We can then view A[ni] as a 2 x 2 matrix with entries ^^^[m]. And again 
we assume that the projection is the same for each physical spin, which is reasonable in practice 
as usual physical systems do have this kind of translational symmetry. The state of the physical 
spins 

\'^) = ([l^j)<S)\^bond)a, (28) 

j a 

where the subscript a labels the bonds, will then have the following form 

|^)= J2 {R\A[zn]---A[zi]\L)\H---Zn). (29) 

Here ij = 0, 1, and \L), \R) are left and write boundary conditions corresponding to the most 
left and most right virtual qubits in Fig.fTOl 

As an example, for the ID cluster state with a projection |0)(00| + which is asso- 

ciated with the bond states \^pcbond), when writing in the form of Eq.([29l). one can equivalently 
choose the projection as 

A[0] = H, and A[l] = HZ. (30) 

The state given by Eq.(l29l) is called a matrix product state (MPS). And it is known that all the 
ID VBS states can be written in this MPS form (for more discussions on MPS, see llSOll '). The 
most famous VBS state is the Affleck-Kennedy-Lieb-Tasaki (AKLT) state, which is the state of 
a spin-1 chain Il2n . The AKLT Hamiltonian of a spin-1 chain is 

haklt ^ }_^^^ . ^ y: i^j:? - 1 (31) 

j j 

Here Sj is the spin operator of the j-th qubit, and -P, j+i is the projection onto the total spin 
J = 2 subspace of each neighboring pair of particles. 

This Hamiltonian is known to be frustration-free and gapped f,21il . If we consider a finite 
chain of n spin-1 particles, then the ground state is four-fold degenerate. One can pick up a 
unique state by appending a spin- 1/2 particle to each end of the chain, and require that the total 
spin of the end spin-1/2 particle and its neighboring spin-1 particle is 1/2. Denote this state by 
\'4^)aklt, then the corresponding Hamiltonian of which \'ip)AKLT is a unique ground state is 

^^i^LT ^ ^ _ ^^^^ _|_ _ Sj^iY + So ■ Si + Sn ■ Sn+1, (32) 

and H^^^'^ is also known to be frustration-free and gapped OTTl . 

The unique ground state of H^^^'^, i.e., the AKLT state \il>)AKLT, has a nice VBS repre- 
sentation [I2TII . as is already mentioned in Sec. HI When choosing the bond state as lipsingiet), 
the projection Pj associated with each circle is to project onto the triplet subspace (e.g the 
symmetric subspace) of two qubits, spanned by 

i=(|00) + |ll)), and -L(|00)-|11)), and -L(|01) + |10)). (33) 

Associated with the bond states \ipbond), the projection matrices based on the MPS representa- 
tion are given by 

A[0] = X, and A[l] = F, and A[2] = Z. (34) 
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As mentioned before, it is already known that the AKLT state \iIj)aklt satisfies Conditions 
2 — 5 in Sec.m Interestingly, it is shown by Brennen and Miyake llITI that \iI))aklt is a re- 
source state for MBQC of a single qubit ll22]| . Their approach is based on the earlier work of 
Gross and Eisert Il20ll . where they developed a framework that shows how an MPS state (and its 
2D generalization) can simulate single-qubit (and two-qubit) gates via single-particle measure- 
ments. This framework allows them to find families of MPS states which satisfy Conditions 
2 — 5 given in Sec.|4]and are resource states for MBQC of a single qubit. These results are very 
nice progress toward finding practical resource states for MBQC, which satisfy all of Condi- 
tions 1 — 5 given in SeclH In this review we adopt another approach that is discussed in Sec.|4l 
i.e., the sufficient condition of reducing the MPS state to a ID cluster state via adaptive local 
measurements at a constant cost. 

In [25], a tabular form of MPS is introduced, which is convenient for showing the reduction 
of other MPS states, including the AKLT state, to a ID cluster state with matrices given by 
Eq.(l30l). In this tabular form, one writes the matrices associated with each physical spin explic- 
itly in a table. Here each column of the table consists of the 2s + 1 matrices of a corresponding 
physical spin of spin s. The physical indices m's determine a selection of the matrices A[m] 
from each column, whose product gives the correct amplitude together with the boundary con- 
ditions {R\ and \L), as given in Eq.(l29l). One then observes from Eq. (l29l) that the following 
properties hold GSlI : 

1. For any two neighboring columns, multiplication of M to the right of all matrices in the 
left column and M^^ to the left of all matrices in the right column simultaneously does 
not change the state. This can be directly seen from the form of Eq. (l29l) . as 

A[i^] ■ ■■A[tk]A[tk-i] ■ ■ ■ = A[t^] ■ ■■A[t,]MAr'A[zk^^] ■ ■ ■ A[z,]. (35) 

2. A unitary transformation in the physical space corresponds to linear combinations of 
entries in the column with coefficients of the unitary. This is a unitary transformation U 
on the basis |m) that results in a new local basis f/|m). 

3. The measurement in the computational basis corresponds to the deletion of column entries 
not consistent with the measurement outcome. This can be directly seen from the form of 
Eq. (|29l). 

4. Columns of a single entry can be removed by absorbing them to a neighboring column. 
This can also be directly seen from the form of Eq. (|29l ). 

As an example. Table 1 of Fig. [T3]is a tabular form that corresponds to two physical spins of 
the AKLT state, which consists of two columns. Starting form this tabular form, we now show 
the reduction from the AKLT state to the ID cluster state, as discussed in [|251. We proceed to 
Table 2 of Fig. [131 which is obtained by adding the Y's with blue colour that represents the same 
state as Table 1 according to property 1 of the tabular form. Table 1 then gives the same state 
as Table 3 of Fig. [131 up to local unitary transformations on the corresponding physical basis, 
according to property 2 of the tabular form. Here in Fig. \T3\ refers to equality up to local 
unitary transformations. As a result, the AKLT state can also be represented by the matrices 
(/, X, Z) as in its MPS representation. 

Starting with this (/, X, Z) form of the AKLT state, we perform two different measure- 
ments Ml and M2 alternatively. Here Mi measures {|0), |i)} versus |2), and M2 measures 
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Figure 13: (This figure is taken from FlG.l in f25\ .) The MPS representation of the AKLT state 
by the matrices (/, X, Z). 



{|0), |2)} versus That is, each measurement consists of a two-dimensional projection and a 
one-dimensional projection. The measurement outcomes are called success (failure) if the out- 
comes correspond to the two(one)-dimensional subspaces, respectively. We measure the two 
measurements sequentially along the AKLT chain, from left to right, and switch the measure- 
ment we use only when the previous one succeeds. 

Table 1 in Fig. [T4l denotes a possible result after these measurements on the AKLT state. 
More specifically, one first measures M.i and succeeds. Next, the measurement M.2 is used. It 
fails and results in the single-dimensional space and we perform it again and it succeeds 
subsequently, etc. After renaming the physical indices and absorbing the X and Z in red color 
to their previous columns, we will get Table 2 in Fig. [14] by properties 4 and 2 of the tabular 
form. This then gives a ID cluster state by the second line of reasoning in Fig. [HI according to 
property 2 of the tabular form. 
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Figure 14: (This figure is taken from FIG. 2 in |[25ll .) The reduction of the AKLT state to the ID 
cluster state. 



The above analysis based on the tabular can be directly generalized to analyze other ID 
resource states. For instance, the modified AKLT state introduced in GUI , which is an MPS 
with 

A[0] = H, and A[l] = X, and A[2] = Y, 

can be similarly shown to be a resource state for MBQC of a single qubit. And the one- 
parameter deformation of the AKLT model considered by Fannes, Nachtergaele and Werner 
in Ref. ffSTl . whose ground state is an MPS with 

A[0] = sm9Z, and = cos ^|0)(1|, and A[2] = cos (0|, 

can also be similarly shown to be a resource state for MBQC of a single qubit fi25l . 
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However, among all these ID resource states for MBQC of a single qubit, the AKLT state 
is of course the most interesting due to its importance in the history of VBS. An experiment 
simulating MBQC on an AKLT state, was performed in an optical system ll52l . In this exper- 
iment, an AKLT state with a single site of spin-1 particle and two end spin- 1/2 particles were 
prepared, and rotations of the initial single-qubit state along any of the X, Y, Z axis via MBQC 
were demonstrated. 

6 A resource state in a spin-5/2 system 

In Sec. [51 we discussed ID resource states which satisfy the natural Conditions 1 — 5 given in 
Sec. m Note that the state lip'y) given by the projection of Eq.(l23l) is a simple example for a 
ID resource state, which is a spin-3/2 state. The advantage of the states discussed in Sec. [51 
including the AKLT state, is that they are spin-1 states. 

To implement universal MBQC, we know that one will need some 2D resource states which 
ideally satisfy the natural Conditions 1 — 5 given in Sec. [H The spin-7/2 state on the honey- 
comb lattice discussed in Sec. [4] provides such an example, but with each physical particle of 
dimension 8. One would wish to reduce the spin dimension, but still keep Conditions 1 — 5 
satisfied. This turns out to be a hard task, as although states satisfying some of the conditions 
may be easy to find, it is in general hard to show that they also satisfy the others. 

The first important step of reducing the dimension is taken by Chen et al. in [[24)1 . They 
constructed a spin-5/2 VBS state on the 2D honeycomb lattice. This 2D honeycomb lattice is 
shown in Fig. [121 where each bond state is chosen as \'ipcbond) as given in Eq.([T7l). Note that a 
2D VBS state is also called a projective entanglement pair state (PEPS) Il53]| . For each spin-5/2 
associated with each circle (with three virtual qubits inside), the corresponding projection onto 
the six-dimensional subspace of the three-qubit Hilbert space is chosen as 

Ptr^c = |0)(000| + |i)(lll| + |2)(100| + |3)(011| + |4)(010| + |5)(101|, (36) 

and the corresponding PEPS state is called the tri-Cluster state, denoted by I'^tric)- 

Based on the discussion in Sec.[4l l^tric) is a resource state for MBQC, as it can be reduced 
to the cluster state associated with the honeycomb lattice via local measurements at a constant 
cost. More precisely, I'^tric) projected onto the subspace spanned by {|0), is the same as 
the cluster state, so are also the states given by {"^tric) projected onto {|2), |3)} and {|4), |5)}, 
up to local Pauli operations. 

Now the task remains is to show that I'^tric) also satisfies Conditions 2 — 5 in Sec.[4l That 
is, one would need to find a Hamiltonian involving only two-body nearest-neighbor interac- 
tions, which is frustration-free, gapped, and has I'^tric) as its unique ground state. We start to 
construct a frustration-free Hamiltonian Htric involving only two-body nearest-neighbor inter- 
actions, which has I'^tric) as its ground state, and then further show that the ground state is 
unique and Htric gapped. 

To construct Htnc, we start from the reduced density matrix of \^tric) for any two neigh- 
boring spins. In general, these two-particle reduced density matrices are dependent on the 
system size, i.e., for different total number of spins n, these two-particle reduced density ma- 
trices are different. However, the range of these two-particle density matrices are independent 
of the system size n. Here the range of a density matrix p is the space spanned by all the 
eigenvectors corresponding to nonzero eigenvalues of p. Therefore, one can find a small system 
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Figure 15: Nearest-neighbor sites of the honeycomb lattice 

to calculate the range of these two-particle density matrices. Because the system has trans- 
lational invariance, we only need to consider three different kinds of neighbors, as shown in 
Figs. [T5al[T5bl[T5cl 

In Fig. [T51 a and b in the circles refer to two types of sites a and b. Note that on the honey- 
comb lattice, the sites of a and b are of different geometry, as one bond on site a goes up and one 
bond on site b goes down. We call the sublattice consisting of all the sites of a type sublattice A 
and the sublattice consisting of all the sites of b type sublattice B. Now denote the two-particle 
reduced density matrices corresponding to the two spins in Figs. ll5alll5blll5cl by pab, Pba and 
Pb respectively, and the corresponding ranges of pab, Pba and pt by Sab, Sba and St respectively. 

a a a 

To compute Sab, as shown in Fig. ll5a[ the virtual qubits 1 to 6 on those sites are only connected 
to virtual qubits a, (3, 7, 5 elsewhere. By tracing out a to 5 from the 5 bonds [tpcbond), we get a 
16-dimensional space for virtual qubits 1 to 6 spanned by \±)i\±)3\^jj singlet) 24\^)5\^)6, where 
|±) = (|0) ±\1))/V2. This 16-dimensional space is then projected by Ptric^ onto qubits 1, 2, 3 
and 4, 5, 6 respectively to give Sab, which is still a 16-dimensional space. Similarly, one can 
compute Sba and St, which are all 16-dimensional spaces. 

Now we choose a two-body Hamiltonian involving only nearest-neighbor interactions as the 
following 

Ht„c = E {Kb + Ka + hi) . (37) 

aGA 

Here the summation is over sites a in sublattice A. And the three terms /i^j,, h^^, correspond, 

a 

respectively, to the projections onto the orthogonal spaces of Sab, Sba and St. Note the overall 
Hilbert space of two spin-5/2 particles is 6 x 6 = 36-dimensional, so h^^, h^^, are all 

a 

projections onto 36 — 16 = 20-dimensional spaces. Apparently, \^tric) is a ground state of 

Htric as K^\^tric) = 0, hlJ-^tric) = 0, and hH'^tric) = 0. Hence iftric is also frustration- 

a 

free. 

To show that I'^tric) is the unique ground state of Htric, we need to verify the condition that 
for any region R of spins in \'^tric), the range Sr of the reduced density matrix on R satisfies 

5r = n Sab ® lR\ab, (38) 

(ab) 

where the intersection is taken over all neighboring pairs ab (i.e., including all the pairs as given 
in Figs. ll5alll5b[[T5cl) . and lR\ab is the full Hilbert space of all spins in region R except a and 
b [[54]| . For every possible configuration containing three or four connected sites in \ '^tric) the 
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condition is confirmed by direct calculation. For larger regions, this condition can be verified by 
induction ll24ll . which is a directly application of Lemma 2 in |[54l . Here we omit the technical 
parts of the results in [[54ll and refer the readers there for more details. 

We now show that Htric is also gapped, following the argument in [l24l. Indeed, there exists 
a constant energy gap t] above the ground state. To estimate the value of t], we first show that r] 
is greater than A, the gap of another Hamiltonian K which also has I'^tric) its unique ground 
state, but involves four-body interactions. 

Consider a Hamiltonian K for a re-labeled version of I'^tric), in which particles are re- 
grouped into disjoint blocks with each containing two nearest neighbors, as shown in Fig. \T6\ 
Note as discussed above, K also has I'^tric) as its unique ground state, where the corresponding 



Figure 16: (This figure is redrawn 



"n"4 • 



Tom FIG. 3 in uM .) Regrouping of lattice sites in tri-Cluster 



state into disjoint blocks, each containing two sites. 



condition similar to Eq. (I38l) can be verified by induction, based on Lemma 2 in [|54]| . 

Now let K = Xlmn ^"1"' where m, n denote two connected blocks, each containing two 
particles m'^l, and n'^l, n^''] respectively. Assuming m'''! and nl'^ are connected, kmn is then 
the projection onto the orthogonal space of the four-particle reduced density matrix of particles 
m''!, m^, nM. We then have 

HtriC = J2^'b>lYl (Kamndrl + ^^MnW + ^nWnw) ^ ^ 5Z Z"^-" = l^"^' ^^^^ 
ab mn mn 

Note that both + ^^w„[;] + h^Wnir]) and k^n are non-negative operators with the 

same null space, so the second inequality holds in Eq. (|39l) for some positive number fi. Without 
loss of generality, one can assume that the gaps of the projectors /i^^ and k^n are both 1, then 
direct calculation gives yU = |. One then has 7] > ^yuA = |A. 

We then bound the gap A. To do this, one can show that K'^ > cK for some positive constant 
c. This is given by direct calculation of and compare it with K, which finally gives c = 1/3. 
Therefore, finally one finds a lower bound on the gap t] of Htric with ?7 > |A > ^. 



7 Resource states in spin-3/2 systems 

After showing that realistic resource states for MBQC could be found in spin-5/2 systems, in 
this section we will see that we could go further, i.e., quantum states of spin-3/2 systems are 
also possible to serve as realistic resource states. 
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7.1 The 2D AKLT state 



In Sec.lH we have shown that the ID AKLT state can be reduced to a ID cluster state, thus it can 
be a resource state of one-qubit MBQC. One would then naturally ask whether any 2D AKLT 
state can be a resource state for MBQC. Recently, Wei et al. and Miyake showed independently 
that the 2D AKLT state on the honeycomb lattice, which is a system composed of spin-3/2 
particles, is a universal resource state for MBQC |[30ll28ll . 

We briefly discuss the approach by Wei et ah. Before starting the discussion, let us explain 
a useful concept, namely quantum encoding (see Chap 10 of ||39l for more details). Remember 
that if a quantum system is two-dimensional, we may call it a qubit. Naturally, if a quantum 
system is composed of more than one qubits, the dimension of its state space will be higher. 
However, it is possible that sometimes only a two-dimensional subspace of the state space of 
this many-qubit system is involved in our consideration. In this situation, for convenience we 
often regard the entire system as one qubit, and we say that the many-qubit system is encoded 
into one encoded qubit or one logical qubit. Similarly, the corresponding Pauli matrices of 
the encoded qubit are called logical Pauli matrices. Furthermore, in a complicated system we 
might encode different subsystems in different ways, and besides, the procedure of encoding 
might also be iterated. As a rule, if a quantum state after encoding is a cluter state, we will call 
this state an encoded cluster state. 

It is discussed in Sec.[5]that the ID AKLT state of spin-l's has a VBS representation. Sim- 
ilarly, 2D AKLT states can also be considered in this manner |[23l . Indeed, 2D AKLT states 
can be defined on any kind of 2D lattices, or even an arbitrary graph [|23]| . For instance, if one 
chooses the projection in Fig. [TT] as the one onto the symmetric subspace of the four virtual 
qubits, the state will be a 2D AKLT state on a square lattice, and the corresponding particles 
are spin-2; if the projection in Fig. [12] is also onto the symmetric subspace, but of the three 
virtual qubits, the state will be a 2D AKLT state on the honeycomb lattice, and the correspond- 
ing particles are spin-3/2. Note that though similar, the two symmetric subspace is different in 
dimension, i.e., the former on the square lattice is of dimension 5 (spin-2), and the latter on the 
honeycomb lattice is of dimension 4 (spin-3/2). 

Let us go back to the work of Wei et al. Based on the observation discussed in Sec. H 
their basic idea is to prove that the 2D AKLT state on the honeycomb lattice can be converted 
to an encoded cluster state on a planar lattice by adaptive local measurements. More con- 
cretely, firstly they found a generalized measurement and perform it on every spin-3/2 particle. 
Secondly, they showed that there exists an encoding scheme determined by the measurement 
outcomes such that the AKLT state on the honeycomb can be regarded an encoded cluster state 
on a random planar graph (i.e., a graph that can be embedded in the plane). Finally, through 
numerical simulation and percolation theory they demonstrated that a typical resulting graph 
state is universal for MBQC. 

As mentioned above and indicated by Fig. [171 in every site (shown as a circle) of the AKLT 
state on the honeycomb lattice C locates one spin-3/2 particle. Being a VBS state, every spin- 
3/2 of £ can be viewed as the four-dimensional symmetric subspace of three virtual qubits. The 
corresponding projection onto this subspace for site v is 

Ps,. = |0)(000| + |i)(lll| + \2){W\ + |3)(iy|, (40) 

where \W) = 7^(|001) + |010) + |100)) and \W) = 7^(|110) + |101) + |011)). Thus, the 
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Figure 17: 2D AKLT state on the honeycomb lattice. Note that all the sites are divided into 
part A and part B. 

AKLT state on £ can be expressed as 

\^AKLt)= Ps,v \i^s^nglet)e, (41) 
v£V(C) e£E{C) 

where V{C) and E{C) are the sets of vertices and edges respectively, and \^singiet)e is a singlet 
state defined in Eq.(fT5l) with its qubits at the two sites connected by e. 
For any site v, consider the following three projections, 

Fv,. = y|(|000)(000| + |lll)(lll|), 

Fv,. = + ++)(+ + +1 + 1 )( 1), 

Fv,y = ^J^{\+i,+i,+i){+i,+i,+i\ + \-i,-i,-i){-i,-i,-i\), (42) 

where |±) = (|0) ± |1))/V2, and \±i) = {\0) ±i\l))/V2. It can be checked that 

Fl^F,^k = Ps,v. (43) 

Namely, these projections form a generalized measurement on each site of the lattice. The key 
part of the work by Wei el al. is that they proved if one performs this generalized measurement 
on every site of the lattice, the resulting state will be an encoded cluster state on some random 
planar graph, up to local unitary operations. Suppose the set of measurement outcomes are 
A = {tty, V E V{C)}, then the resulting state will be a function of A, which is denoted by 
I'^iA)). According to Eq. (|43]), we have 

I^M))= <S) F,,aJ<^AKLT) = (g) F.,a„ (g) l^Psinglet) e- (44) 
v£V{C) veV{C) e£E(C) 
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(a) 




Figure 18: [(b) and (c) of this figure are redrawn from (c) and (d) of FIG. 2 in [|28ll .1 (a) 
Rules Rl and R2. (b) A domain containing two sites, (c) A domain is connected to three other 
domains, and every domain is supported by one site. 

Wei et al. gave a constructive proof to show that |\E'(^)) is an encoded cluster state on some 
random planar graph, up to local unitary operations. Firstly, as a planar graph, the lattice of 
spin-3/2 particles can be recast by two rules illustrated in Fig. fTSl a). These two rules are: (Rl) 
(Edge contraction): Contract all edges e G E{C) that connect sites with the same measurement 
outcomes. (R2) (Mod-2 edge deletion): In the resultant multigraph, delete all edges of even 
multiplicity and convert all edges of odd multiplicity into edges with multiplicity 1 . If several 
sites of £ is contracted into one single site by Rl, we call these sites a domain. Every domain 
will be an encoded qubit in the target encoded cluster state. 

Let us see why a domain can encode a qubit. If the domain contains only one site, the form of 
the measurement suggests that the state of every domain will be in a two-dimensional subspace, 
which naturally means a qubit. The case that one domain contains more than one site can be 
illustrated by an example in Fig.fTSlb). where one domain is composed by two adjacent sites u 
and V with the same measurement outcome a^. Considering the form of F„ ^ in Eq. (I42l) . it can be 
known that the operators Z1Z2, Z2Z2,, Z/^Z^, and Z^Z^ are stabilizer elements of |^(^)). Here 
we use the framework of the stabilizer theory proposed by D. Gottsman |[55ll . For a quantum 
state |a) and a operator A, if A\oi) = \a), we say \a) is stabilized by A, and A is a stabilizer 
element of |«), where A is usually chosen from Pauli operators. Actually, in Sec. [3] we have met 
some examples, and Eq.® in Sec.[3]is one of them, which shows some stabilizer elements of the 
cluster state l^'c). Interestingly, it can be proved that the corresponding cluster state is the only 
state that is stabilized by all the operators of the form given by Eq.® Il39l . which demonstrates 
the purpose of introducing the stabilizer theory, i.e., it is possible to characterize a quantum 
state by its stabilizer elements. Similarly, noting that — Z3Z4 commutes with F„ ^ (g) F^ ^, and 
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that \ip singlet) 'M IS Stabilized by — Z3Z4, one can obtain that — Z3Z4 is also a stabilizer element 
of |\E'(^)). In this way, the state of this domain could be expressed as 

a|(000)„(lll),) + /3|(111)„(000),), (45) 

which therefore might serve as an encoded qubit if one chooses the basis states as |0) = 
|(000)„(111),) and |1) = |(111)„(000),). 

Using the similar idea, we could consider a general domain by the stabilizer theory. Suppose 
there are |C| sites in this domain, and at present we still suppose the corresponding measurement 
outcome is a^. Then |^'(^)) is stabilized by {\i\jZiZj,i, j = 1, 2, 3|C|}. Here = 1 if 
and only if i G v G A, and Aj = — 1 otherwise (note that the 2D AKLT state on the honeycomb 
could be divided into part A and part B as in Fig.[l71). The reason that A/s are introduced is the 
form of the stabilizer elements of \tp singlet), which has a negative sign for ZiZj. As an encoded 

- 3101 

qubit, one can introduce the logical X operation as X = ®j^[Xj, and the logical Z operation 
as Z = XiZi. Similarly, for domains with different measurement outcomes, one can also find 
proper stabilizer elements and logical operations. For example, if the outcome is Uy, one could 
choose the stabilizer generators as {XiXjYiYj, i,j = 1, 2, 3|C|}, X as ^jL[Zj and Z as XiYf, 

if the outcome is a^,., they can be chosen as {XiXjXiXj, i,j = 1, 2, 3|C|}, ®^!fi^j, and AjXj 
respectively. 

When one domain occupies more than one sites, the domain and thus the encoding can be 
simplified. By measurements on those redundant sites, every logical qubit can be supported by 
one site only. For instance, if one measures the site v in the logical qubit of Fig. [TSlb) in the 
basis {1(000)^) ± |(111)„)}, the resulting state will become a |(000)„) 

After determining the encoding scheme, we have to show why the encoded state obtained is 
indeed a cluster state, which can be demonstrated by the example in Fig.fTSlfc). where the state 
is composed by four domains, and each domain contains one site only. Consider the operator 
O = —XiXi/ X2X2' X-^Xy , and we would like to prove that it is in stabilizer of |^'(^)). In 
fact, according to the form of \^jJ)e, it is easy to know that for any i E {1, 2, 3} — XjXj/ is in 
stabilizer of ®eeE{c)H')e- Besides, for any i e {1, 2, 3}, — XjXj/ commutes with for any 
k G {u, V, w}. Therefore, in order to prove is stabilized by C = —XiXitX2X2'XsXy, 

one only need to show that this operator commutes with Fc^a^, which can be verified easily. 

On the other hand, as logical single qubits, we have stipulated logical Pauli operators for 
these four domains, which include Xc = Z1Z2Z3, Z^ = ±-^i', Z^j = ±X2', and Z^, = ±-^3'. 
As a result, it holds that O = iX^Z^Z^Z^. Up to some signs, this is exactly the key character 
of cluster states, as we have introduced in Sec. 3 (see Eq©. 

It should be pointed out that when applying Rule (Rl), multiple edges might be produced. 
As an instance, if there exists a m-multiplicity (m > 1) edge between domains w and v, one 

— — TYi/ — — ( TYi/ mod 2 ) 

will find the resulting state is stabilized by X^Z^ , which is actually X^Z^, because 
of the fact Z^ = I. Therefore, one can find out that the introduction of Rule (R2) will simplify 
the graph while the corresponding resulting state is kept unchanged. 

To prove that 2D AKLT state on the honeycomb is a universal resource state for MBQC after 
getting the encoded cluster states, one still needs to make sure that the connectivity properties 
of the resulting cluster state is good enough. In Ref. Il28]| . Wei et al. provided convincing 
numerical evidences to show that this is indeed the case. 

It should be pointed out that recently Darmawan et al. showed that most states in the dis- 
ordered phase of a large family of Hamiltonians characterized by short-range-correlated VBS 
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Figure 19: (This figure is redrawn from FIG. 1 in ETll .) (a) ID AKLT quasichain. (b) The 
ground state of ID AKLT quasichain is a VBS state. 



states could be reduced to cluster states that are universal for MBQC, and the 2D AKLT state 
on the honeycomb is just an example If56]| . 



7.2 The Quasi-AKLT Chains and Quantum Magnets 

In the subsection above, a nice result that by local operations 2D AKLT states on the honey- 
comb can be converted to 2D graph states is elaborated. As a realistic resource state for MBQC, 
however, one link is still missing for this 2D AKLT. Though widely believed, it remains un- 
known whether the 2D AKLT Hamiltonian on the honeycomb lattice is gapped OTll . That is. 
Condition 4 given in Sec.|4]might not be satisfied. In this subsection, we discuss another work 
by Cai et ah. Based on a new defined system, namely ID AKLT quasichains of mixed spin- 
3/2's and spin-l/2's, they constructed an ideal resource state in a spin-3/2 system for MBQC, 
which successfully satisfies all of Conditions 1 — 5 given in Sec. I4l IISTll . This work was later 
reinterpreted by Wei et al. based on the observation in 4 given in Sec. |4]for universality and the 
technique of generalized measurement discussed in the above subsection ll27l . For convenience 
and consistency, we discuss the approach of ETl to introduce the basic idea of [(57l. 

An AKLT quasichain was defined by Cai et al. as in Fig. [191 a). A little different from 
the original AKLT chain, an AKLT quasichain is composed mainly by some spin-3/2 particles 
forming a backbone, and every spin-3/2 is attached by one or two spin-l/2's. As discussed 
before, every spin-3/2 can be explained as three virtual spin-l/2's followed by the projection 
onto the symmetric subspace, which is shown by Fig. figl b). and the corresponding state can be 
expressed as 

\'^AKLt) ~ ®PSA (g) \^singlet)e, (46) 
A ededge 

where \ipsingiet)e = (|01) — |10))/-y2 is a singlet state on e, and Ps,a is defined as Eq.HOl The 
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Hamiltonian of the quasichain is 



'N~l N 



H=J\Y. ^S... + E + PmI + Pa=1.,. ' (47) 



=1 



where P^h^ is the operator projecting the total spin of a and h onto the spin-/c subspace. It 
can be shown that a AKLT quasichain has a non-degenerate ground state with a finite energy 
gap[|571. 

To be universal for MBQC, 2D model is needed. Based on ID AKLT quasichains, Cai 
et al. constructed an interesting 2D gapped system and proved that its unique ground state is 
universal for MBQC. This 2D model could be illustrated by Fig. |2T1 The key technique in the 
construction is to combine two spin-l/2's from two neighboring ID AKLT quasichains into one 
spin-3/2. The combination could be expressed as a unitary transformation U, 



U= ^ |^,mi + 2m2^^(^i,mi 



mi,m2=itl/2 



1 

bi ' 



(48) 

b2 



which is actually a relabeling of the basis states. Importantly, note that under this unitary trans- 
formation both unitary operations and measurements on spin-l/2's hi and 62 are essentially 
corresponding operations on the spin-3/2 particle B. In this way, the main part of this 2D 
model is spin-3/2 particles, and the only exception is the spin-l/2's on the boundary, which 
actually can also be avoided by considering periodic boundary conditions. The Hamiltonian of 
this 2D model can be expressed as 

H' = ui^H^U\ (49) 

where is the Hamiltonian of the fc-th quasichain, and hi is the tensor product of all necessary 
merging operations given by Eq.(l48l). Cai et al. showed that H' is also gapped [57]. 

By the same generalized measurement given in Eq. (l42l) and similar encoding approach as 
introduced in Sec. l7.1[ the ground state of the AKLT quasichain can be converted to a ID clus- 
ter state by local operations. More precisely, one has to perform the measurement defined by 
Eq.(l42l) on every spin-3/2 particle, then the resulting state turns out to be an encoded cluster 
state, where as discussed before the encoding scheme is determined by the measurement out- 
comes. Meanwhile, a subtle point here is that some spin-l/2's exist. We have mentioned that 
these spin-l/2's are critical in increasing the dimension of the resulting cluster state from 1 to 
2. For this ID case, since the procedure of encoding is almost a repeat of the one in Sec. 17. 1[ 
we provide the encoding details directly. 

As shown in Fig. lSOt a). suppose the measurement outcome on site u is a^, then the logical 
basis states are |0) = |(000)^il2') and |1) = |(lll)u02'), and the logical Pauli operations could 
be chosen as, Z = Zi and X = Xi ® X2 ® X-^ ® X21. Similarly, if the measurement out is ax, 

then these states and operations are |0) = |(+ + +)«(— )2'), |1) = |( )u{+)2'), Z = Xi 

and X = Zx ® Z2 ® Z'i ® Z21 ■ The case of measurement outcome being ay is also similar. 
Note that it is possible that some spin-3/2 is connected to 2 spin-l/2's, as in Fig.[20tb). In this 
case, the encoding is the same except that the state of spin- 1/2 needs to be repeated once. For 
instance, when the measurement outcome is a^, the logical |0) becomes |0) = |(000)tjlol2')- 
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Figure 20: (This figure is redrawn from FIG. 3 in ll27l .) Encoding structure of a logical qubit. 
(a) One spin-3/2 and one spin- 1/2 encode into a logical qubit. (b) One spin-3/2 and two spin- 
1/2's encode into a logical qubit. (c) Two consecutive spin-3/2's and two spin-l/2's attached 
to them encode into a logical qubit. 

The third possibility of encoding a logical qubit is that several consecutive spin-3/2 particles 
have the same measurement outcome, as illustrated in Fig. [20lc). where the two logical basis 
states could be |(000)„l2'(lll)u05/) and |(lll)„02'(000)„l5/). Under the encoding, the ground 
state of every ID AKLT quasichain can be proved to be ID cluster state, which can be verified 
by the stabilizer theory as in Sec. 17.11 Furthermore, since one could simplify every logical qubit 
by measuring the unwanted spin-3/2's and spin-l/2's, it is reasonable to assume that every 
logical qubit occupies only one spin-3/2, which has also be similarly explained in Sec. l7.1[ For 
example, the basis state |0) in Fig. [20la) could be simplified from 1(000)^^12') to |(000)„) by 
measuring qubit 2' in basis (| 0) ± 1 1 ) ) / V^. 

After knowing how to obtain logical ID cluster states, the main challenge remaining is to 
produce a 2D cluster state by somehow connecting them together. In order to overcome this, 
Cai et al. proposed the following interesting solution, as illustrated by Fig.[2Tl 

To produce universal resource states for MBQC, for arbitrary two neighboring logical qubits 
u and V on different quasichains, one needs the capability to perform two possible operations 
between them: One is a logical controUed-Z gate defined in Eq.®, and the other is to keep the 
state of every logical qubit unchanged. Cai et al. proved that both of them could be realized by 
operations on qubits hi and 62 only, up to local operations on logical qubits u and v. Recall that 
unitary operators and measurements on qubits hi and 62 are essentially corresponding operations 
on the spin-3/2 particle B. 

Without loss of generality, we assume that when producing ID cluster states, both of the 
measurement outcomes on sites u and v are a^. To keep the two logical qubits unchanged, one 
first applies two X operations on qubits hi and 62. then measures them in basis (|0) ± \\))/\/2 
separately. If the new measurement outcomes are mi and m2, it can be checked that the effective 
operation is to perform and Z"'"^ on logical qubits u and v respectively. Actually, this is 
to essentially simplify the encoding of logical qubits u and v by removing the unwanted spin- 
1/2's, and the resulting state might need the correction of local unitary operations. 

Let us turn to the case of controUed-Z operation. Once more, in this case two X operations 
are needed on qubits 61 and 62 firstly, then a controUed-Z gate is applied on these two qubits. 
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Figure 21 : 2D model based on ID AKLT quasichains. Spin- 1/2 pairs in the circles are mapped 
into spin-3/2's. 

Up to local operations, a logical controUed-Z gate is implemented. 

In short, we see that ID AKLT quasichains can be connected together and then be converted 
to a logical 2D cluster state, which is universal for MBQC. By showing that the operations 
needed for related spin-l/2's correspond to operations on the resulting spin-3/2's state through 
the mapping in Eq.(l48l). we see that that the ground state of the 2D model defined by Cai et al. 
can be a realistic resource state for MBQC. 



In the above sections, we have seen successful examples showing the existence of spin-5/2 
and spin-3/2 systems that are ideal resource states for MBQC. As the most practical quantum 
systems, it is natural to ask whether it is possible to find ideal resource states for MBQC from 
qubit systems. Unfortunately, it has been pointed out by Chen et al. that this is not the case (32^. 

The basic idea of Chen et al. is to show that for any two-body frustration-free qubit Hamil- 
tonian, there always exists a ground state which is a product of single- or two-qubit states. That 
is to say, any entangled state cannot be the non-degenerate ground state of such a Hamiltonian, 
which rules out the possibility of qubit systems being ideal resource states for MBQC. 

Consider an n-qubit state |^) which is genuinely entangled. That is, |^) cannot be a product 
state with respect to any bipartition of these n qubits. We denote the reduced density matrix of 
qubits i and j by pij. Consider the following Hamiltonian if^ 



where Uij is the projection onto the kernel of pij. 

It can be observed that is a two-body frustration-free qubit Hamiltonian, and |\[') is 
in the ground space of iJ^. Actually, one can find that S'(^) is the smallest among 

all the ground-state spaces of two-body frustration-free qubit Hamiltonians that contain \'^). 
For a general two-body frustration-free qubit Hamiltonian H = J2k ^k, Hk's might not be 
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projections. However, we could replace them with n'^,s and the ground-state space does not 
change, where Uk is the projection onto the orthogonal subspace of the ground-state space of 
Hk. In this way, we only need to consider Hamiltonians composed of projections, which are of 
the form of Eq.dSO]). 

Next, let us note an important point that helps with the proof. Suppose there exist 2x2 
nonsingular linear operators Li, L„, and |\&) = where L = Li (g) ■ ■ ■ ® Ln. In the 

language of quantum information theory, two states |$) and |^') can be transformed to each 
other via L are called equivalent under stochastic local operation and classical communication 
(SLOCC)[[58l. Note that |$) is a ground state of H if and only if |$) is a ground state of 

H' = ^(L, ® L,)tn,,(L, ® L,), (51) 

and that if one of |\1') and |$) is product state, the other one also is. Thus, if one could show 
there exists a product state in the ground-state space of two-body frustration-free Hamiltonian 
containing |^), the similar situation will happen for any two-body frustration-free qubit Hamil- 
tonian having |$) as a ground state. Therefore, we only need to consider the SLOCC equivalent 
states with respect to the transformation L. 

We start from the simple case of three qubits. It is known that there only two different 
SLOCC equivalent classes for three-qubit genuinely entangled states [[591 , which are repre- 
sented by the W state 

|W^) = -^(|001) + |010) + |100)), (52) 

and the GHZ state 

|Gi/Z) = i=(|000) + |lll)). (53) 
It is straightforward to check that 

5(|W^)) = span{|iy),|000)}; (54) 

and 

SdGi/Z)) = span{|000), (55) 

Apparently, both S'(|Vr)) and S{\GHZ)) contain at least one product state. 

Next we proceed to the four-qubit case. Since |\1') is a genuinely entangled state, the rank 
of every Iljj is at most 2. We then discuss two cases: all pij are of rank 3 or 4, or at least one pij 
is of rank 2. 

If all Pij are of rank 3 or 4, then is actually the so-called homogenous Hamiltonian 
discussed in [[60][ . According to Lemma 2 of [[60][ . one can always find a product of single-qubit 
state in the ground- state space of H^. We refer the readers to [[601 for technical details on the 
discussion of homogenous Hamiltonians. 

If for some pair of qubits, the rank of pij is 2. Without loss of generality, we assume {i, j} 
to be {3, 4}. Now we can reduce the problem to a three-qubit case by encoding qubits 3 and 4 as 
one single qubit. Suppose P34 is supported on two orthogonal states 1^0)34 and 1^1)34. Consider 
the following isometry, 

^: 10)3' ^|^/'o)34, 11)3' ^1^^1)34, (56) 
where qubits 3 and 4 are encoded by a new qubit 3'. 
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Define 

\^) = V^\^),H = V^H^V. (57) 

It can be seen that H is also a two-body frustration-free Hamiltonian, and |$) is still genuinely 
entangled. Besides, |\E') is a ground state of ifi^ if and only if |$) is a ground state of H. 

Since H is exactly the case of three qubits we have already discussed, it follows immediately 
that the ground-state space of i/' contains a product state |ai)®|a2)®|a3'). Lety|a3') = 1(3^)- 
If 1/^34) is product state, then |^) is also a product state. For the case that I/334) is entangled, one 
can always find a product state ® |/32) in the range of P34, since this is a two-dimensional 
subspace of two qubits ll6T]| . 

The case ofn>A qubits can just be similarly shown by induction. Thus we finally conclude 
that for any entangled state |^), there always exists a product state of single qubits which is 
also a ground state of Hq,. This also indicates that if we do not require that |^) being genuinely 
entangled, then the ground-state space of i?^ always contains a state which is a product of 
single- or two-qubit states. 

Therefore, we can conclude that any genuinely entangled n-qubit state cannot be a 
unique ground state of a two-body frustration-free Hamiltonian. Because that a resource state 
for MBQC must be genuinely entangled. Conditions 1 — 5 given in Sec.|4]cannot be simultane- 
ously satisfied for any qubit system. In other words, there does not exist any idea resource state 
for MBQC in spin-1/2 systems. 

Input: a set of product 
states that spans the ground 
space of the simplified ho- 
mogenous Hamiltonian. 

Figure 22: (FIG.l in (331.) The general structure of the ground space of 2-body frustration-free 
qubit Hamiltonian 

It should be pointed out that the above discussion regarding the ground-state space proper- 
ties of two-body frustration-free Hamiltonians of qubits can be strengthened further. In fact, it 
is shown by Ji et al. that the structure of the ground-state space of this kind of Hamiltonians 
can be characterized completely [[33ll . That is, such a ground-state space is a span of tree tensor 
network states of the same tree structure, which is illustrated in Fig. [22] (we refer the readers 
to Il48ll for detailed discussion tree tensor network states). Here we briefly discuss the structure 
given by Fig. [22] as follows. For any two-body frustration-free Hamiltonian H of qubits, one 
can produce a corresponding special Hamiltonian H', called the simplified homogenous Hamil- 
tonian (again, we refer the readers to [[6OI for technical details on the discussion of simplified 
homogenous Hamiltonians). It can be shown that the ground-state spaces of these simplified 
homogenous Hamiltonians are actually spanned by (not necessarily orthogonal) product states 
of single qubits. In Fig. [221 every blue triangle stands for an isometry defined in Eq. [56] The 
very left side of the tree structure is the input, and it is the set of the product states that span 
the ground-state space of the corresponding simplified homogenous Hamiltonian H'. Then ev- 
ery product state solution goes through the tree structure from the left to the right, and after 
being operated by the meeting blue triangles, each of them becomes a solution of the original 
Hamiltonian H. Finally, the entire ground-state space of H is then spanned by these tree tensor 
network states. 
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9 Future directions 



We have reviewed the recent resuhs toward finding practical resource states for MBQC. Start- 
ing from Verstraete and Cirac's formulation in 2004 [|26ll relating cluster states to VBS states, 
we discuss a general method for producing VBS resource states for MBQC by Gross and Eis- 
ert [|20l . There has been good progress made to identify practical resource states for MBQC 
with VBS states. One-dimensional resource states that are useful for single-qubit MBQC are 
identified for spin-1 systems, as well as two-dimensional resource states for spin-5/2 and then 
later spin-3/2 systems. 

Down the road there are still many interesting open threads. The no-go result by Chen et 
al. discussed in Sec. [Drules out the existence of any spin-1/2 system as a practical resource 
state for realizing MBQC, satisfying all of Conditions 1 — 5 given in Sec. HI As the resource 
states satisfying all of Conditions 1 — 5 are identified in Sec.|7]only for spin-3/2 systems, this 
naturally leaves a gap on whether a resource state satisfying all of Conditions 1 — 5 can be 
identified in a spin-1 system. In other words, it is still unknown whether there exists a spin-1 
state on some two-dimensional lattice which could be universal for MBQC, and at the same time 
is the unique ground state of some frustration-free and gapped Hamiltonian involving only two- 
body nearest-neighbor interactions. So far such a state has only been found for one-dimensional 
lattices, which does not rule out the existence of such a state on a two-dimensional lattice. 

Given that the spin-1 /2 systems are the most commonly found systems in nature, one would 
also like to seek for some relaxation of the ideal Conditions 1 — 5 for the systems. For instance, 
one may try to relax the "frustration-free" condition (Condition 5) or the "uniqueness" condition 
(Condition 4) requirement by resorting to some other physical mechanisms, such as topological 
protection as discussed in Il34l[35]| or perturbation as discussed in Il36ll37l[38l . 

In a more general theoretical framework, each one of Conditions 1 — 5 lacks a deeper un- 
derstanding. That is why the construction of each resource state discussed in this review turns 
out be somewhat "aJ hoc", apart from the fact that they are all VBS states. 

For Condition 1, we would like a general understanding concerning which kinds of states 
could be universal for MBQC. In particular, whether the observation discussed in Sec. |4l that 
is, if a spin state I?/;) can be reduced to a resource cluster state \'^c) via adaptive local mea- 
surements at a constant cost, then l^/^) is a resource state for MBQC, is in general a necessary 
condition for all resource states. So far, all the examples of resource states discussed in this 
review satisfy this condition, and indeed we use this condition to prove the universality prop- 
erty for MBQC for all these states. It then remains open whether a state which can be used for 
universal MBQC, but does not satisfy this condition, can be identified. 

For Conditions 2 and 3, we would like a general understanding that which kinds of quantum 
states can be unique ground states of two-body Hamiltonians of certain interaction patterns 
(for instance, two-body nearest-neighbor interactions on some kind of lattices). Recently, this 
problem has been linked to the certain kind of correlations of quantum states by Chen et al. 
in Il62l . The work along the line of Chen et al. raises a general question of the so-called "from 
ground states to local Hamiltonians". That is, to develop general understanding and method to 
find some desired local Hamiltonians to have a given quantum state as its unique ground state. 

However, these local Hamiltonians constructed by Chen et al. in 1621 are in general frus- 
trated, and their discussion is restricted to finite systems (i.e., systems composed of finite n 
particles), so whether these Hamiltonians are gapped, which is Condition 4, needs further dis- 
cussion. It remains a general challenge to determine whether a two-body Hamiltonian is gapped. 
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and one would like to obtain concrete examples showing whether the 2D AKLT state on the 
honeycomb lattice is gapped, as mentioned in Sec. l7.1[ 

Finally, we would like to mention that Condition 5, i.e., the study of ground states of 
frustration-free Hamiltonians, is closely related to computer science. Indeed, some of the tech- 
niques used in Sec. [H have been used to deal with some quantum version of a computer science 
problem in II6OII . This problem, called the quantum 2-satisfiability (2-SAT) problem, is of vital 
importance in the study of the theory of quantum computational complexity. Indeed, the charac- 
terization of the general structure of the ground-state spaces of two-body frustration-free qubit 
Hamiltonians shown by Fig. [22] has an immediate and interesting corollary regarding quantum 
computational complexity theory, i.e., the corresponding counting problem of quantum 2-SAT 
is equal to its classical counterpart in computational complexity, which answers an open prob- 
lem raised in ^63^ . More details can be found in If33ll . 

We therefore believe that these recently progresses related to MBQC with VBS states dis- 
cussed in this review will not only push the effort of realizing large scale quantum computers, 
but also open up new directions that can further enhance the research in both quantum informa- 
tion science and condensed matter physics. 
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